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Abstract: In order to predict and analyse the potential geographical distributions of three important quarantine 
invasive pests, Dacus bivittatus, D. ciliatus and D. vertebratus, three ecological niche modeling techniques, 
BIOCLIM, DOMAIN and Maximum entropy niche-based modeling ( Maxent ) were implemented by using 
distribution records of the three fruit fly species and a set of environmental predictor variables. Differences in 
prediction performance of the three models with thresholds were observed. An evaluation using independent 
records of D. bivittatus showed that Maxent offers the most accurate predictions than two other models based on 
three values of ROC/AUC, Kappa, and TSS. Prediction outcomes made by Maxent revealed that the three fruit fly 
species have broadly similar potential ranges in Central American, South American, Southeast Asia, and Coastal 
areas of Australia in general. D. ciliatus has the comparatively widest potential range among the three species, 
including Coastal areas of Mediterranean Sea, Saudi Arabia, Yemen, Oman and South Iran, suggesting that it may 
be tolerant of the widest range of climatic conditions among the three species. In China, while large areas of 
Yunnan and Hainan are very habitable for all the three fruit fly species, southern part of Guangdong and Taiwan 
are also their habitable areas. D. ciliatus has the widest potential distribution area, with southern part of Sichuan, 
Guizhou and Tibet plus the coastal areas of southern China all being its suitable areas. The risk of the three fruit 
fly species permanently establishing in Guangdong if introduced exist but low. Jackknife analysis revealed that 
temperature and its variation have comparative significant influence on the distribution patterns of three fruit fly 
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species both in global and restricted regions. 
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1 INTRODUCTION 


Guangdong province is situated at southern China 
where the climate is suitable for fruit fly species, and is 
adjacent to Hong Kong. It is a major gateway for 
international commercial trade and traveler 
transportation in China mainland. About one third of 
total imported fruit consignments of the country enter 
China mainland through Guangdong ports each year. 
From 2003 to 2007, Dacus bivittatus ( Bigot) and D. 
ciliatus Loew had been intercepted four and three times 
respectively in traveler inspection at Guangdong ports, 
and there are other five times of Dacus spp. interception 
records in addition (unpublished data) (Liang et al., 
2006, 2007). The Chinese Government classified Dacus 
spp. as the important quarantine invasive pests because 
of the economic and international business significance 


(AQSIQ, 2005; CABI, 2005). D. bivittatus, D. 


ciliatus and D. vertebratus Bezz are three important and 


serious invasive pest species of Dacus spp., whose 
native distribution areas are Africa ( Liang et al., 1993; 
CABI, 2005; Liang et al., 2006; EPPO, 2007), and 
there are no distribution reports of them in China yet. 
Nowadays, with more and more frequent business trade 
between China and African countries, the interception 
frequency and potential risk of introduction of Dacus 
spp. into China are growing higher. Given the tropical 
and subtropical climatic features of southern China, and 
the wide areas of fruit, vegetable and forest there, it is 
very likely that the Dacus spp. pests cause serious 
damage to the crops and environmental security if they 
introduced ( Liang et al., 2006, 2007). Researches 
about predictions of potential distributions of D. 
bivittatus, D. ciliatus and D. vertebrates are few but 
necessary. 

Therefore, objectives of this study were two: (1) 
Predicting and comparing potential geographical 
distributions of D. bivittatus, D. ciliatus and D. 
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vertebratus, both in China and global area. The process 
was launched based on ecological niche requirements 
(largely in climatic dimensions) and their ecological 
distributions in their origin regions. Models of potential 
distribution are developed and tested quantitatively for 
D. bivittatus as a measure of the predictive power of 
the methodology, then applied to two others; (2) 
Finding the critical environmental variables determining 
the distribution patterns of the three fruit fly species. 


2 MATERIALS AND METHODS 


2.1 Occurrence data 

Native-range distributional data for D. bivittatus , 
D. ciliatus and D. 
Royal Museum for Central Africa ( http :// projects. 
bebif. be/enbi/fruitfly ). The 
Afrotropical specimen data for the ceratitidine 
( Tephritidae, Dacinae, tribe Ceratitidini ) fruit fly 
genera Ceratitis, Trirhithrum Bezzi, and Capparimyia 


vertebratus were obtained from 


database contains 


Bezz, and summarizes the holdings of all major 
European, African and North American natural history 
collections, supplemented with a few records from 
specimens housed in the National Insect Collection at 
Pretoria, South Africa ( De Meyer et al., 2008). 
Non-native distributional data was gathered from 
literature, research papers, reports, and electronic 
resources. The data reflect both contemporary as well 
as historical records of the three species outside Africa, 
and permanent as well as temporary establishments. 
Only records with authoritative confirmation of 
those with doubtful 
identification were excluded from model building, 


identification were included; 


besides we used only records of permanent 
establishment for distribution predicting. In total, 345, 
343 and 243 records were adopted for D. bivittatus , 
D. ciliatus and D. vertebratus , respectively. 

Localities were georeferenced based either on 
coordinates provided directly on specimen labels for 
individual specimens, or on coordinates extracted from 
official gazetteers ( e. g., GeoNet, http://gnswww. 
Records with 


unspecified or unknown localities were removed from 


nga. mil/geonames/GNS/index. jsp). 
consideration. In cases of multiple records for a 
particular locality, only one was used in analysis. 
2. 2 Environmental data sets and relevant 
geographic layers 

The environmental layers “ bioclimatic’ were used 
in predicting procedure, which consist of important 
ecological factors and climatic features worldwide. It 
contains 19 layers with 2. 5?spatial resolution, and was 
derived by the WorldClim project, http ://biogeo. 
berkeley. edu. These were Annual Mean Temperature, 
Annual Precipitation, Isothermality, Max Temperature 
of Warmest Month, Mean Diurnal Range, 
Temperature of Coldest Quarter, Mean Temperature of 


Mean 


Driest Quarter, of Warmest 


Quarter, Mean Temperature of Wettest Quarter, Min 


Mean ‘Temperature 


Temperature of Coldest Month, Precipitation of Coldest 
Quarter, Precipitation of Driest Month, Precipitation of 
Driest Quarter, Precipitation of Warmest Quarter, 
Precipitation of Wettest Month, Precipitation of Wettest 


Quarter, Precipitation Seasonality, Temperature 
Annual Range and Temperature Seasonality. 
Procedures producing these interpolated climate 


surfaces are described by Hijmans et al. (2005). 

Layers of the 1:4 million-scale about administrative 
boundary of China were downloaded from National 
Fundamental Geographic Information System ( NFGIS: 
http ://nfgis. nsdi. gov. cn/). The layers characterized 
the geographic location of provinces boundary, main 
rivers, important cities, large habitat areas, main road 
lines and main railway lines of China. We covered them 
on the final prediction outputs to analyse the distribution 
patterns of three fruit fly species. 

No vegetation-cover layers and architectural 
composition layers were used, because they change 
markedly year by year; and besides given the 
heterogeneous nature of habitats ( including man-made 
horticultural environments ) could potentially be 
occupied by the target fruit fly species because of their 
polyphagous host range ( De Meyer et al., 2008). 

2.3 Ecological niche modeling 

The comparison of prediction performance were 
made among three methods, i. e., BIOCLIM ( Busby, 
1991), DOMAIN ( Carpenter et al., 1993) and 
Maxent ( Phillips et al., 2006). The final prediction 
maps of three fruit fly species were produced based on 
species occurrence data using the superior one. In 
each method, some rules or mathematical algorithms 
were used to define the ecological niche of the species 
based on the distribution records of the fruit fly species 
in the multidimensional environmental space. Once the 
species niche is defined, its projection into the 
geographical space produces a predictive map. 

BIOCLIM defines the ecological niche of a species 
as the bounding hyper-box that encloses all the records 
of the species in the climatic space (Busby, 1991). 
in the 


environmental space, defined by the most extreme 


Thus, it creates a rectilinear ' envelope ’ 
(minimum and maximum) records of the species on 
each environmental variable. To reduce the sensitivity 
of model predictions to outliers, the species records are 
sorted along each variable and only the records that lie 
within a certain percentile range of the data are used 
for model construction. 

DOMAIN uses a point-to-point similarity metric 
(based on the Gower distance) to assign a value of 
habitat suitability to each potential site based on its 
proximity in the environmental space to the closest 
(most similar) occurrence location ( Carpenter et al., 
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1993). A threshold value of suitability can then be 
selected to determine the boundaries of the ecological 
niche. In contrast to all previous methods, 
environmental envelopes defined by DOMAIN are not 
necessarily continuous in the environmental space. 
BIOCLIM and DOMAIN were launched on DIVA-GIS 
Version 5.4.0.1 ( Hijmans et al., 2001) , which was 
downloaded from http://www. diva-gis. org/. We 
followed recent research papers to tune the parameters 
and thresholds of the two models, and used the True- 
False option, to place pixels into one of two categories. 

In this study, for BIOCLIM, two percentile ranges 
were applied: 90% (disregarding 5% of the values on 
each side) , 95% (disregarding 2.5% of the values on 
each side ); for DOMAIN, two thresholds were 
chosen: 0. 95 and 0. 97 (Hernandez et al., 2006; 
Beaumont et al., 2007; Tsoar et al., 2007). 

Maxent is a machine learning method, which 
estimates the distribution of a species by finding the 
probability distribution of maximum entropy (i.e. that is 
closest to uniform) subject to constraints representing our 
incomplete information about the distribution. The 
constraints are the expected value of each environmental 
variable should match its average over sampling locations 
derived from environmental layers. The model evaluates 
the suitability of each grid cell as a function of 
environmental variables at that cell. The suitability values 
provided by Maxent range from 0 (unsuitable habitat) to 
100 ( optimal habitat ). Maxent version 3. 2. 19 
downloaded from http://www. cs. princeton. edu/ ~ 
schapire/maxent/ was used to do the predictions. We 
followed recently published best practice approaches to 
tune its parameters (Phillips et al., 2006; Phillips and 
Dudik, 2008). The default settings were used for most 
Random Test 
percentage: 0, Regularization multiplier: 1, Convergence 
threshold: 0. 00001, Logistic output format) , and with a 
few exceptions 1 000). In 
addition, Maxent has a built-in jackknife option through 


Maxent parameters ( Auto Features, 


( maximum iterations : 


which the importance of individual environmental data 
layers can be estimated (Fig. 5). 
2. 4 Accuracy assessment and performance 
comparison methods 

Maxent predictions ( Maximum training sensitivity 
plus specificity logistic threshold ) were compared with 
two thresholds of predictions of BIOCLIM and DOMAIN 
respectively in this study. 
( Mann-Whitney U) was used to determine whether one 


A non-parametric test 


prediction is significantly better than another in the 
same sample sets. Three methods were used to assess 
the prediction accuracy of models: First, the area 
under the curve ( AUC) in receiver operating 
characteristic analyses (ROC) (Elith, 2002; Phillips 
et al., 2006); Second, the Cohen’ s kappa values 
from an error matrix (Cohen, 1960) ; Third, the True 
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Skill Statistic (TSS) ( Allouche et al., 2006). All 
layer processing was done on ArcGIS 9. 2 Desktop 
(ESRI Inc., 2006). 

Receiver operating characteristic (ROC) plot was 
selected as the first method to evaluate the accuracy of 
generated distributions. A curve that maximizes 
sensitivity for low values of the false-positive fraction is 
considered a good prediction and is quantified by 
calculating the AUC. The AUC can be used as a 
measure of overall performance of the model and has 
values usually ranging from 0.5 (random) to 1. 0 
( perfect discrimination ) , if the values below this range 
indicating a prediction that is worse than random 
(Engler et al., 2004). The main advantage of ROC 
analysis is the area under the ROC curve ( AUC) 
providing a single measure of model performance, 
independent of any particular choice of threshold. 
While generally used when data are available, ROC 
plots can also be created with presence and background 
absence data. In this study, Hawth' s Analysis Tools 
for ESRI’ s ArcGIS (v3. x) software ( Beyer, 2007) 
was implemented to randomly generate 100 background 
observations in global; and we put them into the ROC 
curve procedure in place of absence points to 
accompany the presence observations ( Hernandez et 
al., 2006). The AUC values are calculated in SPSS 
16.0 (SPSS Inc., 2007). 

Cohen’ s kappa is the most popular measure for the 
accuracy of presence-absence predictions ( Loiselle et 
al., 2003; Petit et al., 2003; Rouget et al., 2004; 
Segurado and Araújo, 2004). The kappa statistic ranges 
from —1 to +1, where +1 indicates perfect agreement 
and values of zero or less indicate a performance no 
better than random. Both commission and omission 
errors are accounted for in one parameter. 

Although both ROC/ AUC and Cohen’ s kappa are 
widely used in prediction accuracy assessment, several 
studies have criticized their associated shortcomings 
( Allouche et al., 2006; Lobo et al., 2008). So true 
skill statistic (TSS) was adopted as the alternative 
measure of prediction accuracy assessment which is not 
affected by prevalence and the size of the validation 
set, and two methods of equal performance have equal 
TSS scores ( Allouche et al., 2006). Like kappa, TSS 
takes into account both omission and commission 
errors, and success as a result of random guessing, 
and ranges from — 1 to +1, where + 1 indicates 
perfect agreement and values of zero or less indicate a 
performance no better than random. 

2.5  Jackknife test 

A jackknife test was then performed for each 
species to assess the relative contribution. of the 
individual environmental variables in the prediction 
procedure. The test was performed as follows: (i) the 
full complement of 19 environmental variables was used 
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to perform the prediction as a baseline for comparisons ; 
(ii) each environmental variable was eliminated 
systematically to provide the 19 possible combinations 
of 18 environmental variables each, and the 19 
possible combinations were used to analyze the 
importance of each missing variable; (iii) similarly, 
but in the opposite way, each environmental variable 
was included systematically in single-variable analysis 
(Peterson and Cohoon, 1999; Levine et al., 2007). 

Although jackknife tests can be performed by 
hand, we used Maxent to conduct them because the 
whole procedure was automatically executed, and the 
results were intuitive and easy to interpret ( Phillips, 
2006; Phillips et al., 2006). 
2. 6 Application procedures for predictive 
performance comparison and potential 
distribution prediction 

The whole approach comprised of six main steps 
(Fig. 1): (1) Dividing the occurrence samples of D. 
bivittatus into two parts stochastically: 80% and 20% 
parts by a SAS program (SAS Institute Inc., 2003 ) ; 
(2) Using the 80% of D. bivittatus occurrence points 
to create the prediction maps by BIOCLIM, DOMAIN 
and Maxent; (3) Testing the prediction accuracy of 
three models by ROC/AUC, Cohen' s kappa and TSS; 
(4) Repeating Step 1 to Step 3 ten times; (5) 
Finding out the superior model; (6) Using the all 
occurrence samples of each fruit fly species to predict 
their habitat areas by the superior model. Since our 
goal is to predict global potential habitat areas, we test 
model performance only for D. bivittatus (8096 of the 
records) , which has the most occurrence records ( De 
Meyer et al., 2008). In total, each partition for D. 
bivittatus held 276 training points and 69 test points. 
Ten random partitions were made to assess the average 
behavior of models. 


3 RESULTS 


3.1 Performance comparison 

Fig. 2 indicates that Maxent has the best 
performance (with the highest AUC, Kappa and TSS) ; 
DOMAIN has the median performance; and BIOCLIM 
with both thresholds performed comparatively poorly. 
DOMAIN with 0.97 threshold has the closest AUC value 
with Maxent, and the difference is not significant (two- 
tailed Mann-Whitney U test: P of AUC = 0. 65); 
DOMAIN with 0.95 threshold has the closest Kappa and 
TSS value with Maxent, and the differences are both 
Extremely significant (two-tailed Mann-Whitney U test: 
P of Kappa = 0.000, P of TSS = 0.000). Based on 
the performance comparisons mentioned above, we used 
the predictions of Maxent to analyze the potential 
distributions of the three species of fruit flies, which 
present each pixel’ s probability in the form of 
continuous value ( corresponding logistic value). 
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Fig. 1 Procedures for the predictive performance 
comparison among BIOCLIM, DOMAIN and Maxent 
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Models & Thresholds 
Fig. 2 AUC, Kappa and TSS value comparisons among 
BIOCLIM, DOMAIN and Maxent in different thresholds 
The error bars represent 95% confidence interval for mean. On X axis, 
the figures beside * BIOCLIM' and ' DOMAIN’ are their percentile 


ranges, the threshold of Maxent is the maximum training sensitivity plus 


Maxent 


specificity logistic threshold. 
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3.2 Prediction of potential distribution areas 

Fig. 3 shows the potential distribution areas of 
three fruit fly species in globalarea made by Maxent. 
Although the difference of global area distribution 
coverage of three fruit fly species is not significant, D. 
ciliatus has the comparatively widest distribution areas 
among three fruit fly species, the latter is D. 
bivittatus; and D. vertebratus has the 
distribution areas. For D. ciliatus, the long island of 
USA, wide areas of Mexico, Bahamas, Cuba, Central 
Africa, Saudi Arabia, 


Yemen, Oman, south Iran and Sri Lanka, some areas 


smallest 


America, South America, 


of Europe (costal areas of Portugal, Spain, Italy and 
Greece), Yemen, Jordan, Oman, Pakistan, India, 
Southeast Asia and Australia (coastal areas) are all its 


potentially suitable habitat areas ( Fig. 3: B). In 
addition , it is the only species among three species that 
has the potential distribution in coastal areas around 
the Mediterranean Sea and Southern Europe. For D. 
bivittatus , the long island of USA, Mexico, Bahamas, 
Cuba, Central America, middle of South America, 
Africa, Yemen, Sr Lanka, some areas of Pakistan, 
India, Southeast Asia and Australia ( coastal areas) 
are its suitable habitat areas ( Fig. 3: A). For D. 
vertebratus, the long island of USA, wide areas of 
Mexico, Bahamas, Cuba, Central America, South 
America, sub-Saharan Africa, Yemen, Sri Lanka, 
some areas of Pakistan, India, Southeast Asia and 
Australia (coastal areas) are its suitable habitat areas 
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Fig. 3 Global potentially suitable areas of three fruit fly species predicted by Maxent and presented each 
pixel' s probability in the form of continuous value ( Corresponding logistic value) 
A: D. bivittatus; B: D. ciliatus; C: D. vertebratus. The darker color of a pixel indicates the higher confidence in fruit fly presence. White dots on maps 


represent the locations and distributions of the three fruit fly species. 
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Fig. 4 Potential distributions of three fruit fly species in China predicted by Maxent and presented each pixel’ s probability in 


the form of continuous value ( Corresponding logistic value ) 
A: D. bivittatus; B: D. ciliatus; C: D. vertebratus. The darker color of a pixel indicates the higher confidence in fruit fly presence, and bold lines on maps 
represent the Yangtze River. 
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Fig. 4 and Fig. 5 reveal that large areas of 
Yunnan and Hainan provinces of China are very 
habitable for all three fruit fly species, southern part of 
Guangdong and Taiwan are their habitat areas too. It is 
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indicated that D. ciliatus has the widest potential 


distribution areas in China, including southern part of 
Sichuan, Guizhou and Tibet. 
southern China are its suitable areas. 


The coastal areas of 








008 £f 


Fig. 5 Potential habitat areas of three fruit fly species in Guangdong predicted by Maxent and presented each pixel’ s 
probability in the form of continuous value ( Corresponding logistic value) 
A: D. bivittatus; B: D. ciliatus; C: D. vertebratus. The darker color of a pixel indicated the higher confidence in fruit fly presence. The areas covered 
by 1096 simple hatch with 45? angle is Guangdong province, covered by 1096 simple hatch with 90? angle is Fujian province, and covered by 1096 
crosshatch is Hainan province. White dots on the map are important cities of Guangdong where entry-exit ports are located ( AQSIQ, 2005). 


3.3  Jackknife analysis 

Fig. 6 indicates that results of jackknife tests of 
three fruit fly species are similar with each others on 
the whole, and five variables have notable influences 
on distribution prediction. These are Isothermality, 
Mean Temperature of Coldest Quarter, Min 
Temperature of Coldest Month, Temperature Annual 
Range, and Temperature Seasonality. Meanwhile, five 
variables show less contribution to the prediction since 
they achieve few gains. These are Mean Diurnal 
Range, Precipitation Seasonality , 
Coldest Quarter, Precipitation of Driest Month, and 


Precipitation of 


a 


Annual Mean Temperature 

Annual Precipitation 

Isothermality 

Max Temperature of Warmest Month 
Mean Diurnal Range 

Mean Temperature of Coldest Quarter 
Mean Temperature of Driest Quarter 
Mean Temperature of Warmest Quarter 
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Precipitation of Driest Month 
Precipitation of Driest Quarter 
Precipitation of Warmest Quarter 
Precipitation of Wettest Month 
Precipitation of Wettest Quarter 
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E With all variables 








= 
2.0 0.00.2 0.40.60.8 1.01.2 1.41.61.80.0 0.5 1.0 1.5 2.0 
Regularized Training Gain 


Precipitation of Driest Quarter. 

For D. bivittatus, D. ciliatus, and D. vertebratus 
in that order, the environmental variables with highest 
gain when used in isolation are Temperature 
Seasonality , Isothermality and Isothermality, 
respectively (the longest blue bars in Fig. 6: A, 
B, C). 

In addition, the light blue bars in Fig. 6 reveal 
that for all the three kinds of fruit fly, no variable has 
more useful information that is not already contained in 
the others, because omitting each one in turn didn’ t 
decrease the training gain much. 

b C 





Regulanzed Training Gain Regularized Training Gain 


Fig. 6 Jackknife tests for environmental variable significance performed by Maxent 
A; D. bivittatus; B: D. ciliatus; C: D. vertebratus. 
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4 DISCUSSION 


Performance comparison of ecological niche 
models in this study indicates that Maxent has the best 
performance among the three models (Fig. 2) ; and in 
former similar studies, Maxent turned out to be the 
best when compared with other methods ( ENFA, 
GARP, Maxent and PCA), which makes it a very 
promising method ( Segurado and Araújo, 2004; 
Graham and Hijmans, 2006; Tsoar et al., 2007). 
Besides, there are two additional advantages make it 
outstanding. First, Maxent can be easily interpreted by 
property 
importance. Its construction idea is to estimate the 
target distribution by finding the distribution of 


maximum entropy (i. e. that is closest to uniform ) 


human experts, a of great practical 


subject to the constraint that the expected value of each 
feature under this estimated distribution matches its 
empirical average. This turns out to be equivalent, 
under convex duality, to finding the maximum 
likelihood Gibbs distribution. For species distribution 
modeling, the occurrence localities of the species serve 
as the sample points, the geographical region of 
interest is the space on which this distribution is 
defined, and the features are the environmental 
variables ( Phillips et al., 2004). Second, it has a 
high success rates and statistical significance with 
sample sizes as low as five occurrence records ( Pearson 
et al., 2007). 

The prediction accuracy of GARP and Maxent is 
greatly influenced by exactly which observations are 
included; the importance of each record depends on 
whether that locality represents a unique environment 
not represented by the other sample points ( Pearson et 
al., 2007). In this research, the distribution records 
of each fruit fly species covered all its native-range 
distribution countries and included the important non- 
native distribution areas, which represent a wide 
variety of environment conditions and make the final 
outcomes be of high accuracy. 

The final prediction outcomes indicate that the 
risk of three fruit fly species permanently established in 
Guangdong exists but low (Figs. 4 and5). D. ciliatus 
has the widest potential distribution areas in China, 
followed by D. bivittatus and D. vertebratus. Wide 
ciliatus. In 
addition, the neighbor countries of China (e. g., Laos, 
Thailand and Vietnam) , where are also suitable for the 


areas of Fujian are suitable for D. 


three fruit fly species, should not be ignored. 

Based on jackknife analysis, temperature and its 
variation influence the distribution patterns of three 
frut fly species significantly both in global and 
restricted regions ( Fig. 6), and which is consistent 
with former research ( Vayssieres et al., 2008 ). 


Temperature is a very important abiotic factor affecting 
immature development of Dacus spp. Knowledge of 
relationships between developmental time — and 
temperature is essential for predicting population 
fluctuations in the field. According to Vayssieres 
(2008), there was a strong linear regression model 
between temperature and development rate of D. 
ciliatus. Temperature thresholds (T, ^C) and thermal 
constants (K, degree days) of the different immature 
ciliatus ( egg, larva, pupa and total 
developmental stage in order) on pumpkin are: 
10.56°C and 26 degree days, 10. 23*C and 90 degree 
days, 11. 20°C and 207 degree days, and 11. 04°C 
and 323 degree days, respectively. 


Since 1980, the African people have emerged into 


stages of D. 


China ( Guangdong mostly) at a high speed which 
keeps raising intensely per year. In 2005, the number 
of incoming travelers from Africa at Guangzhou is up to 
31 766, and the average increasing rate per year is 
37.996. 
African in Guangzhou and Shenzhen ranges from 1/50 
to 1/25 ( Li et al., 2008; Pan et al., 2008). It is very 
possible that D. bivittatus, D. ciliatus and D. 
vertebratus China with 


Nowadays, the population probability of 


introduced into incoming 
travelers. 

Many biological invasions result from human 
activity. They are mainly considered as one of the 
consequences of increasing international trade. Once 
an exotic fruit fly species is established, it cannot, in 
most cases, be eradicated. Attempts at management 
are often not particularly effective ( Vayssieres et al., 
2008 ). 
current issue for both tropical and temperate countries. 
In all especially 
temperature, are likely to affect the invasion process. 


Forecasting tephritid invasion is a major 


cases, environmental factors, 
Also, these results are useful for improving laboratory- 
rearing methods and for building simulation models to 
predict fruit fly population dynamics. 

The process of invasive species entering and 
establishing in a new region is complex, requiring 
about six barriers according to Richardson (2004) . A, 
main geographical barriers; B, local environmental 
dispersal 


( disturbed 


habitats); F, environmental barriers ( undisturbed 


barriers; C, reproductive barriers; D, 


barriers; E, environmental barriers 
habitats ). In this study, potential habitat areas in 
China ( barrier B) for the three species of fruit flies 
were estimated. We should also investigate potential 
competition of the invasive exotic species with 
indigenous species or with other introduced species, 
which is considered an important part of barriers C and 
D (Duyck et al., 2004; De Meyer et al., 2008 ). 
Today Bactrocera dorsalis ( Hendel) and Bactrocera 
cucurbitae ( Diptera: Tephritidae) are two most serious 


fruit fly species in China, and caused enormous 
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damage to fruit industry every year. They are widely 
distributed in southern provinces of China ( Hou and 
Zhang, 2005; Zhou et al., 2007). Further relevant 
studies should concentrated on the competition between 
invasive species and indigenous species based on 
concept of propagule pressure and Allee effects ( Drake 
et al., 2005; Lockwood et al., 2005; Drake and 
Lodge, 2006). 
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基于 最 大 燃 值 法 生态 位 模型 (Maxent) 的 


三 种 实 蝇 潜在 适 
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(1. 中 山大 学 有 害 生 物 控制 与 资源 利用 国家 重点 实验 室 /昆虫 学 研究 所 , J M 510275; 
2. 广东 出 和 人 境 检验 检疫 局 检验 检疫 技术 中 心 , 广州 510623) 


摘要 : 本 研究 首先 对 3 种 重要 生态 位 模型 BIOCLIM, DOMAIN 和 Maxent ( JT sc 48 (E JE BR) 的 分 布 预测 精确 度 
进行 了 分 析 和 比较 , 再 结合 分 布点 记录 以 及 一 系列 环境 数据 图 层 对 3 种 重要 外 来 人 侵 性 检疫 害虫 ( 萌 芦 赛 紧 实 晶 
Dacus bivittatus IRAE WW 3E SEC D. ciliatus VS JI XE RK D. vertebratus) 的 潜在 适 生 性 分 布 区 域 进行 了 预测 和 
分 析 。 在 模型 预测 精确 度 的 比较 过 程 中 , 3 BEDE TRES ROC/AUC, Kappa, TSS ) 均 显示 Maxent 拥有 最 好 的 预测 结果 


和 最 好 的 运行 


性 能 。 由 Maxent XY ij XE 38 RAE AR E ERE SCBRCNTPS JR EE ESSE Wa H3 BRUM AR. o, 这 3 种 实 蝇 


在 中 美洲 、 南 美洲 、 东 南亚 和 澳大利亚 沿岸 的 广大 地 区 在 总 体 上 具有 相似 的 分 布 区 域 。 相 对 而 言 , REAR SEE 
蝇 在 全 球 范围 具有 最 为 广泛 的 分 布 区 域 , 除 前 述 地 区 外 ,其 潜在 适 生 区 还 包括 地 中 海 沿 岸 、 沙 特 阿拉 伯 、 也 门 、 安 曼 


和 伊朗 南部 的 大 片 地 区 , 这 也 意味 着 在 3 RRR, 


它 能 忍受 变化 幅度 最 广 的 生态 .环境 条 件 。 在 中 国 ， 云 南 和 


海南 都 极 适宜 于 3 种 实 蝇 的 生存 , 同时 广东 南部 及 台湾 的 部 分 地 区 也 是 它们 的 潜在 适 生 区 。 基 于 Maent 的 预测 结 
果 显 示 , 相对 而 言 , REESE SEP ERA AR AT IZ Kk, 除 前 述 省 份 和 地 区 外 ， 四川 贵州 
和 西藏 的 南部 部 分 地 区 以 及 中 国 南 部 的 部 分 沿海 地 区 , 也 都 是 它 的 潜在 适 生 区 。 综 合 所 得 出 的 预测 结果 , 3 ER 
实 蝇 从 境外 传人 广东 并 在 此 和 定 殖 的 风险 可 能 性 是 实际 存在 的 。Jackknife 分 析 显 示 , 温度 以 及 与 此 有 关 的 环境 因子 对 
于 3 种 实 蝇 在 全 球 和 局 部 地 区 的 分 布 模式 和 分 布 情况 都 有 极 大 的 影响 , 并 需要 进一步 的 研究 。 
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